Motility states of molecular motors engaged in a stochastic 
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Abstract Intracellular transport is mediated by molecular motors that pull cargos along cy- 
toskeletal filaments. Many cargos move bidirectionally and are transported by two teams 
of motors which move into opposite directions along the filament. We have recently intro- 
duced a stochastic tug-of-war model for this situation. This model describes the motion of 
the cargo as a Markov process on a two-dimensional state space defined by the numbers of 
active plus and active minus motors. In spite of its simplicity, this tug-of-war model leads 
to a complex dependence of the cargo motility on the motor parameters. We present new 
numerical results for the dependence on the number of involved motors. In addition, we de- 
rive a simple and intuitive sharp maxima approximation, from which one obtains the cargo 
motility state from only four simple inequalities. This approach provides a fast and reliable 
method to determine the cargo motility for a given experimental system. 

Keywords Molecular motors ■ intracellular traffic ■ bidirectional movement • stochastic 
processes ■ nonequilibrium 

1 Introduction 

Molecular motors are cellular proteins which transform the free energy released from chem- 
ical reactions into mechanical work. Cytoskeletal motors walk in a directed fashion along 
cytoskeletal filaments and thereby transport cargo through the cell [1]. Examples are kinesin 
and dynein, which walk along microtubule filaments, and certain myosins, which walk along 
actin filaments. Single motors of this kind have been studied extensively in recent years both 
experimentally [39] and theoretically [23]. The filaments possess an intrinsic direction: they 
have one 'plus' and one 'minus' end. They are essentially one-way roads for the motors. For 
example, cytoplasmic dynein walks to the microtubule minus end, while kinesin- 1 walks to 
the microtubule plus end. Most cells have a unidirectional microtubule cytoskeleton [20]: 
The microtubule minus ends are typically located near the cell center, while the plus ends 
point outwards towards the cell periphery. In special cells such as neurons with their long 
axonal protrusions, the microtubules form a unipolar parallel array with the minus ends 
pointing towards the cell center and the plus ends pointing towards the axon tip. 

Although each motor walks only into one direction along such an isopolar filament 
network, many cellular cargos, like mitochondria, pigment granules and endosomes, are ob- 
served to move bidirectionally, reversing direction every few seconds [42]. Therefore, both 
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plus-end and minus-end motors must be involved in the transport of a single cargo, and in- 
deed both kinesin and dynein are found simultaneously on various cellular cargos [12]. We 
have recently proposed a stochastic tug-of-war model for this situation, in which the mo- 
tors are coupled via the mechanical interaction with their common cargo [29]. This model 
maps the motion of the cargo to a stochastic process on a finite two-dimensional state space 
defined by the numbers of active plus and active minus motors. Despite the simplicity of 
the model, the cargo exhibits a rather complex motility behavior. Depending on the motor 
parameters, the cargo motion exhibits several qualitatively different motility states, which 
correspond to qualitatively distinct steady state solutions of the corresponding Markov pro- 
cess. In particular, for biologically relevant parameter values, the cargo exhibits switching 
between fast plus and fast minus motion with or without interspersed pauses, as found in 
experiments [12,42]. This fast bidirectional motion, which is usually associated with a co- 
ordination mechanism rather than with a tug-of-war, is obtained in our model via a dynamic 
instability, which leads to a high probability of having only one motor type bound at a given 
time. Therefore, the tug-of-war is a cooperative mechanism for bidirectional transport. In 
addition to the fast bidirectional motion, our tug-of-war model can explain many properties 
of bidirectional cellular cargo transport observed in experiments in a quantitative way [29] . 
This implies that the signalling pathways that control intracellular transport may directly 
target the different motor molecules rather than an additional coordination complex [29,43]. 

In our previous study in [29], we have emphasized the experimental relevance of our 
model and given a detailed comparison between theory and experiment. In this article, we 
describe the modeling procedure in detail and present new numerical results for the param- 
eter dependence of the motility states. The dependence on the motor numbers leads to the 
prediction that bidirectional cargo transport with pauses can only be accomplished by at 
least 4 motors on the cargo. Furthermore, we describe a simple and intuitive approximation, 
which we call the 'sharp maxima approximation' and which provides a reliable description 
for the dependence of the different motility states on the single motor parameters. Using the 
analytical solution obtained for this approximation one can easily obtain the prediction of 
our model for the motion of a cargo earned by two motor species in a simple and transparent 
manner without performing tedious numerical calculations. The approximation also allows 
to calculate the effect of a change in the motor properties, which might be accomplished 
e.g. by changing the concentration of ATP or salt, by adding regulatory molecules such as 
dynactin, or by mutation of the motors. 

Cooperative effects in systems of many molecular motors have been studied extensively 
in the theoretical literature [2,5,7, 11, 15, 16, 18,24,41]. Typically, these studies have con- 
sidered large numbers of motors belonging to a single motor species as appropriate, e.g. for 
the modeling of muscles. In contrast, transport in cells, as mentioned, is often bidirectional 
because of the cooperation of (at least) two motor species and is typically based on rather 
small numbers of motors, in the range of 1-10 [14, 18]. Cargo transport over cellular length 
scales, both unidirectional and bidirectional, has been modeled in a coarse-grained way as 
movement of a particle with effective rates characterizing the speeds and run lengths of the 
cargo [23, 24, 27, 31, 35]. While these models provide an appropriate description of cargo 
movements on large time and length scales, they cannot address the stochastic fluctuations 
of motor-filament binding and unbinding, the effect of the motors on each other, and the de- 
pendence of cargo transport on the molecular properties of the motors, which are addressed 
within the framework of our model. 

This article is organized as follows. We introduce our stochastic tug-of-war model in 
Section|2]and discuss its numerical solutions in Section[3] These two sections give a detailed 
description of the model introduced in [29], and summarize the main theoretical results 
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Fig. 1 Cargo with a total number of N + = 2 plus (blue) and N- = 2 minus (yellow) motors, (a) The numbers 
«+ and n_ of plus and minus motors that actually pull the cargo, which are (n+,re_) = (2, 1) in the figure, 
fluctuate because the motors stochastically bind to and unbind from the filament, (b) Two-dimensional state 
space of (N+ + l)(N- + 1) = 9 states labeled by (« + ,n_). 

reported in [29]. In addition, new numerical results for varying motor numbers are shown. 
In Section |4] we present a new analytic 'sharp maxima' approximation in order to explain 
this parameter dependence. This approximation allows to determine the motility states of a 
given system by using only four simple inequalities. 



2 Modeling 

We consider a cargo particle with a fixed number of N + plus and AL minus motors. Typically 
these numbers will be in the range of 1 to 10 motors as observed for many cargos in vivo 
[14, 18]. For N + = or N- = 0, we recover the model for cooperative transport by a single 
motor species as studied in [18]. Each of the motors can bind to, move along, and unbind 
from the filament. Because of the stochastic binding and unbinding, at each time only n + 
plus and «_ minus motors are bound to the filament, with < n + < N + and < «_ < N-, 
see Fig.[TJa). Since only these bound motors can exert force on the cargo, the cargo motion 
is determined by the numbers (n + ,n_). Unbinding or binding of a plus motor in the state 
(n_|_,«_) occurs with rate e + (« + ,«_) or 7T + (n + ,«_) and changes this state to (« + — l,n_) 
or (« + + respectively. Analogously, unbinding or binding of a minus motors occurs 

with rate £_(«+,«_) or tt_(w + ,«_), respectively, This leads to a 'random walk' on the state 
space shown in Fig. |Ttb). This random walk or Markov process is described by a Master 
equation for the probability p(n + ,n-,t) to have n + active plus and «_ active minus motors 
at time t . This Master equation has the form 

^-p(n + ,n-,t) = p(n+ + l,«-,f) £+(«+ + l,n_) + p(n + ,n- + 1, ?)£_(«+,«_ + 1) 
at 

+ p(n + — \,tl-,t) lt+{n+ — l,n_) + p(n + ,«_ — l,t) 7T_(ra + ,n_ — 1) 
—p(n + ,n-,t) \n + (n + ,n-) + 7T_(n + ,«_) +£ + (n+,n_) +£_(«+, n_)] .(1) 

In order to obtain the rates £+(_) and 7T+(_) for unbinding and binding, respectively, of 
one motor, we first establish a model for a single motor, which is based on the transport 
properties of molecular motors as measured in single-molecule experiments. We then derive 
the effective cargo rates that enter the Master equation |Q~|| by assuming that the motors feel 
each other only because they act on their common cargo, as explained below. 
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Parameter 


symbol 


kinesin-1 


cytoplasmic dynein 


stall force 


F s 


6pN [34,36] 


7 pN ('strong') [37] 


detachment force 


F d 


3 pN [34] 


1.1 pN ('weak') [26,44] 
? 


unbinding rate 


£o 


1/s [34,40] 


0.25/s [17,32] 


binding rate 


»0 


5/s [3,21] 


1.5/s [6,40] 


forward velocity 


Vf 


1 nm/s [6,40] 


1 jllm/s [17,37] 


superstall velocity amplitude 


VB 


6 nm/s [6] 


? 



Table 1 Single motor parameters of our model, and values for kinesin- 1 and cytoplasmic dynein, taken from 
the cited references. A question mark indicates that this parameter is not available. 



2.1 Model for a single motor 

A single motor can bind to a filament with rate K, walk along it with velocity v, and unbind 
from it with rate e. In our motor tug-of-war, opposing motors exert force on each other, so 
that the load-dependence of these parameters is crucial. For some molecular motors such 
as kinesin-1, this load-dependence has been measured in single molecule experiments. The 
motor velocity v decreases with increasing load force from its zero-load value Vf to zero at 
the so-called stall force F s . This decrease is approximately linear, and given by 

v(F) = vf (1 -F/F,) for < F < F„ (2) 

as has been observed for kinesin-1 [6,30,36], for kinesin- 3 [38] and cytoplasmic dynein [37]. 
Here we use the convention that the load force F is positive if it acts opposite to the motor's 
forward direction. For superstall load forces F > F s , the motor can walk backwards, but 
only slowly, as has been shown for kinesin-1 [6,30] and cytoplasmic dynein [19,25]. In this 
regime, the functional form of the motor force-velocity relation is unclear; for simplicity, 
we use a linear form 

v(F) = v B (l-F/F s ) for F>F S . (3) 

The superstall velocity amplitude vb characterizes the slope vb /F s of the force-velocity re- 
lation in the superstall regime. Our results are essentially independent of the exact form of 
the motor force-velocity-curve as long as it decreases monotonically with external force and 
exhibits slow backward motion. The unbinding e rate of the motor, which equals £o at zero 
load force, increases exponentially with the load force as 

e(F) = e exp[F/F i ], (4) 

which defines the detachment force Fj. Such a functional form has been measured for 
kinesin-1 [36], and also follows theoretically from Kramers or Bell theory [4]. The force- 
dependence of the binding rate is difficult to assess experimentally. However, it is expected 
to depend only weakly on the load force because an unbound motor relaxes and then binds 
from its relaxed state. We therefore take the binding rate equal to 

7i(F) = Kq (5) 

for all load forces F . 

All single-motor parameters have been measured for the motor kinesin- 1 , see Table [T] 
Other motors have parameters of similar orders of magnitude, but their parameters are not as 
well established. For cytoplasmic dynein, for example, no measurements for the detachment 
force and for the superstall velocity amplitude are available, and different labs have reported 
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different values for the stall force, see Table[T] In this work, we use the kinesin-1 parameters 
for the missing dynein parameters, and we call dyneins with stall force 7 pN 'strong' and 
with 1.1 pN 'weak' dyneins. 

2.2 Effective rates for the cargo 

We assume that the motors on the cargo act independently and feel each other only because 
(i) opposing motors act as load, and (ii) same-directional motors share this load. If each plus 
motor feels the load F + (and generates the force — F + ) and each minus motor feels the load 
— F_ (and generates the force F_), this means that the force balance on a cargo pulled by n + 
plus and «_ minus motors is 

n+F + = -n-F-=F c (n+,n-), (6) 

where F c denotes the force which acts on the team of the plus as well as on the team of the 
minus motors. Here, the sign of the force is chosen positive if it is a load on the plus motors, 
i.e. if it points into the minus direction. If only one motor type is bound, i.e. if n+ = or 
«_ = 0, then F + = F_ = F c = 0. The force balance as given by © represents Newton's third 
law: each motor feels the same force as it exerts. A single bound plus motor thus feels the 
force F + = F c /n + . Using this in the single motor unbinding rate Q, this implies that the 
effective rate for the unbinding of one plus motor in the state (ra + ,w_) is 

£+(«+,«_) =«+£<)+ sxp[F c (n+,n-)/(n + F d+ )]. (7) 

Here and below, the index '+' labels plus motor properties, e.g. F d+ is the detachment force 
of a single plus motor. The prefactor n + in equation l|7} describes that there are n + bound 
plus motors available for unbinding. Analogously, the effective rate for the binding of one 
plus motor in the state where there are N + — n + unbound plus motors, is, with 

equation (O, 

7f+ (n+, nJ) = (N+ - n+ ) Uq+ . (8) 

Analogous expressions hold for the unbinding and binding rate of a minus motor with the 
parameters indexed by '— '. 

The cargo force is determined by the condition that both plus motors, which experience 
the force F+, and minus motors, which experience the force — F-, move with the same 
velocity v c as given by 

v c (« + ,«_)=v + (F + ) = -v_(-F_). (9) 

Here, the sign of the velocity is taken to be positive in the plus direction and negative in 
the minus direction. The velocity balance relation l|9} assumes that all motors walk with the 
same velocity. This can be justified by considering what happens when this steady state is 
disturbed, e.g. by the binding of a minus motor while the cargo is steadily moving into the 
plus direction. The 'new' motor has bound in a relaxed state and feels a low load force. It 
therefore steps forward, gets stretched and burdens itself with part of the load force imposed 
by the plus motors, taking away part of the load force of the other minus motors. It does 
so until all minus motors feel roughly the same load force. This justifies our mean field 
treatment that each plus or minus motor feels on average the same force F + or — F_ . The 
new force balance depends on whether the plus or minus motors are 'stronger'. If the plus 
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motors together can produce a larger force than the minus motors, i.e. if n + F s+ > w_F,_ 
(with the 'new' number n_ of bound minus motors), then the new minus motor runs forward 
until it feels a superstall load force and then starts to run backwards. The total force on the 
cargo increases, and all motors slow down until a new force balance is achieved. If the minus 
motors can produce a larger force than the plus motors, i.e. if h__F,_ > n + F s+ , the minus 
motors 'take over': the new minus motor runs forward and takes away more and more load 
force from the other minus motors until they all feel a substall load force — F_ < F s _ . The 
plus motors, on the contrary, start to move backward as soon as they feel a superstall load 
force F + > F s+ . The forward and backward stepping of motors thus provides a relaxation 
mechanism which establishes the force and velocity balances as provided by equations $6^ 
and (9). 

The force and velocity balances as given by equation $6$ and lead to the cargo force 

F c (n + ,n_) = X{n+,n-)n + F s+ + [\ -A(n+,n_)] «_F S _, (10) 

with 

A(n+,«_) = 1/ [1 + (n + F s+ vo-)/(n^v +)} , (11) 
and to the cargo velocity 

I \ n+F s+ -n-F s ^ 

n+F s+ /v + +n_F s _/vo- 

Here the single motor velocity parameter vo+ is equal to the plus motor forward velocity vf+ 
or superstall velocity amplitude vb+, depending on whether the plus motors move forward, 
i.e. v c > 0, or backward, i.e. v c < 0, respectively. Analogously, the single motor velocity 
parameter vo- is equal to vf- for v c < and to vb- for v c > 0. It follows from equation i l 12b 
that the cargo direction is determined by the intuitive 'majority rule' 

v c > for n + F s+ > n_F s _, i.e. plus motors 'win' 

v c < for n + F s+ < n__F s _, i.e. minus motors 'win' (13) 

v c = for n + F s+ = n_F s _, i.e. 'tie'. 

Since < A < 1, the cargo force as in equation d 1 Ob is a convex combination of the maximal 
plus and minus motor stall forces n + F s+ and it-F s -. Therefore, if for example the plus 
motors win, these motors walk under substall and the minus motors under superstall load 
force, i.e. ti-F s - < F c (n + ,n_) < n + F s+ , which is consistent with the single motor force- 
velocity relation equations ® and (3]t. 

The cargo force given by equation d 1 Ob determines the plus motor rates ^ and ([8} and 
the corresponding minus motor rates for the Master equation 1Q. However, one has to take 
care of reflecting boundary conditions, which ensure that the numbers (n + ,n_) of active 
motors stay within the intervals < n + < N + and < n_ < AL. 

Experiments in vivo usually monitor cargos which have been walking along a filament 
for some (unknown) time. This has two implications for our calculations: First, we are in- 
terested in long-time properties of cargo motion, i.e. in the time-independent steady state 
probability p(n + ,n_) to have n + active plus and n_ active minus motors, which is obtained 
by setting the time-derivative in the Master equation (Q equal to zero. Second, we are inter- 
ested in probabilities that are conditioned on the cargo being bound to the filament. We can 
separate cargo unbinding into two steps: (i) the last motor unbinds from the filament so that 
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the cargo is in the state (w + ,w_) = (0,0) but still close to the filament, and (ii) the cargo dif- 
fuses away from the filament and then undergoes free diffusion in the surrounding solution. 
We call this diffusing state U. The rates for cargo unbinding into and rebinding from the 
state U depend on the geometry of the system and the viscosity of the surrounding solution. 
In Appendix lAl we show that the probability which solves the time-independent 

Master equation (TJ is in fact the steady-state probability conditioned on the cargo being 
bound to the filament, i.e. not to be in state U, and that this probability is independent of 
the rates that connect the state (0,0) with state U. We can therefore ignore cargo unbind- 
ing and diffusion. Furthermore, a cargo in the state (0, 0) is close to the filament. Binding 
of plus and minus motors from this state is then described by the intrinsic motor rates 7Zq + 
and 7io_, respectively, which justifies our choice of the binding rates in $Q also for the state 
(«+,«_) = (0,0). 

We solve the Master equation (Q for the steady state by determining the eigenvector 
of the associated transition matrix with eigenvalue zero. In addition, we simulate individual 
cargo trajectories by using the Gillespie algorithm [10] for the binding/unbinding dynamics 
as given by equations ^ and l[8} and let the cargo move with velocity v c in the intervals 
between (un-)binding events. 

The Master equation (Q} describes a two-dimensional random walk on the network 
shown in Fig. Gib). The steady-state solution of the Master equation (Q} is a nonequilib- 
rium steady state, as can be seen by considering the ratio of the product of forward and the 
product of backward rates of the cycles of the network of Fig.[TJb). This product is not equal 
to 1 for all cycles, as would be required for an equilibrium state [22,23]. The rates of the 
Master equation are nonlinear in the state space variables (n + ,rc_) because of the nonlinear 
force-dependence of the unbinding rates (J7J, which leads to a cooperative effect, namely 
the unbinding cascade described in the next section. The emergence of cooperative behavior 
arising from the nonlinear force dependence of the unbinding rate has also been proposed as 
an explanation for collective effects in muscles [7, 15] and mitotic spindle oscillations [1 1]. 



3 Motility states 

Depending on the number and the parameters of the motors on the cargo, the cargo ex- 
hibits qualitatively different types of motions which we call 'motility states'. We define the 
motility states via the number and locations of the maxima of the motor number probability 
p(n + , «_ ) in the steady state. These maxima characterize the cargo motion because the cargo 
spends most of its time in these maxima (n + ,n_). With respect to cargo motion, there are 
three different types of maxima: a maximum at (n + ,n_) with n + > 0, n_ = corresponds 
to fast plus motion with cargo velocity v c = Vf+, a maximum at («+,«_) with n + = 0, 
ra_ > corresponds to fast minus motion with cargo velocity v c = v/r_, and a maximum 
at («+,«_) with n + , n_ > corresponds to small cargo velocity v c . The latter can be seen 
from equation d 1 2b by considering the fact that biological motors have a small superstall 
velocity amplitude Vb <S vf- Expanding the cargo velocity v e (« + ,«_) in, e.g. for winning 
plus motors with n + F s+ > ri-F s -, vb-/\>f+, leads to v c « \>b- \n+F s+ / (n-F s J) — 1], which 
is small for small vb- Intuitively, the losing minus motors have to walk backwards, which 
they do only slowly, so that cargo motion is slowed down correspondingly. 
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Fig. 2 Possible realizations of the symmetric tug-of-war. (a) The tetrameric motor Eg5, walking along an- 
tiparallel microtubules in the mitotic spindle, corresponds to a symmetric tug-of-war of N^- = 1 'plus' motor 
(the upper half of the motor) and N- = 1 'minus' motor (the lower half of the motor), (b) A large cargo carried 
by four motors of the same type along two antiparallel microtubules corresponds to a symmetric tug-of-war 
of N + = 2 'plus' (upper) and N- = 2 'minus' (lower) motors. 

3.1 Symmetric tug-of-war 

It is instructive to consider the 'symmetric' tug-of-war for the same number N = N + = N- 
of plus and minus motors with identical single motor parameters apart from their forward 
directions. This case is theoretically appealing because of its simplicity and intuitive appeal. 
In addition, the effects of the motor tug-of-war are most pronounced in the symmetric case 
because the motors are 'equally strong'. Moreover, this symmetric situation may not be too 
unrealistic for transport in vivo where plus end and minus end transport are often found to 
exhibit astonishingly similar properties [13,31,44]. Furthermore, the most relevant plus and 
minus motors, kinesin-1 and cytoplasmic dynein, exhibit similar properties in terms of their 
processivities [17,34] and their force- velocity-curves [37], although the latter is still contro- 
versial, see e.g. [9]. The symmetric tug-of-war can also be realized by using only one type 
of motor but antiparallel filaments, see Fig. [2] The tetrameric kinesin-5 motor Eg5 walks 
along antiparallel microtubules in the mitotic spindle during cell division. This corresponds 
to a symmetric tug-of-war of N+ = 1 'plus' motor, corresponding to one dimer, and N- = 1 
'minus' motor, corresponding to the other dimer of the Eg5 tetramer, see Fig.|2ja). In a pos- 
sible in vitro experiment, a large cargo is transported by several motors of one type along 
antiparallel microtubules, see Fig.[2jb). This again corresponds to a symmetric tug-of-war of 
'plus' and 'minus' motors, where the assignment of a motor to the 'plus' or 'minus' motor 
type depends on the microtubule to which it attaches. 

In the symmetric situation, the indices for plus and minus motors for the single motor 
parameters can be omitted, i. e.: 

7R) = = Xq-,£o = £o+ = Bo- , F s = F s+ = F s - , F d = F d+ = F d _ , v = v 0+ = v _ (14) 

In this symmetric case, inspection of the Master equation equation 10 and its rates ([7]l and 
<[8j> shows that the steady state probability p(n+,nJ) only depends on four dimensionless 
parameters: 

the motor number N = N+ = the desorption constant K = Eq/hq, 

the force ratio / = F s /F d , and the velocity ratio v = vb/vf. (15) 

As most motors walk backwards rather slowly, the velocity ratio is very small for most 
motors, e. g. v = 0.006 for kinesin-1. Because of the small value of v, the results are rather 
insensitive to its precise value and very close to the results for no backward motion, i. e. 
v = 0. For fixed numbers N of motors on the cargo, the relevant parameters are therefore 
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Fig. 3 The motility states for the symmetric tug-of-war of N = N + = N- = 5 plus and minus motors are 
characterized by qualitatively different motor number probability distributions (top row) and cargo trajectories 
(bottom row), as described in the text. Motor parameters are for kinesin as given in Table \T\ except for (a) 
F s = 0AF d , £q = 0.2^b, (b) £q = 0.35^b, (c) F s = l.QFj, £q = Q.2hq. In all cases the superstall velocity 
amplitude is scaled as vg = lnm/s-F s /pN. 

the force ratio / and the desorption constant K. The latter parameter describes the binding 
affinity of the motors; for processive motors like kinesin or dynein it is smaller than 1. If 
the force ratio / ' = F s /F^ > 1, the force F s that the motor exerts exceeds the force Fj that 
it can sustain. This ratio is most important for the qualitative properties of the steady state 
solution p(n+,nJ) of the Master equation (Q. As described in [29], there are three types of 
solutions, which we call motility states, and which differ in the number and locations of the 
maxima of p(n + ,n_), see Fig.[3]and Fig. [4] 

(0) No motion: For 'weak motors' with small stall to detachment force ratios / = F s /Fd, 
the motors are hardly affected by the presence of the opposing motors and bind / unbind as 
if they were alone. The probability distribution /?(«+,«_) has a single symmetric maximum 
at a state with n + = see Fig.[3]al), corresponding to a state with zero velocity. The cargo 
therefore shows only small fluctuations around its starting position, see Fig.[3ja2). 

( — h) Fast bidirectional motion: For 'strong motors' with large force ratio /, the mo- 
tors feel the opposing motors strongly, and spontaneous symmetry breaking occurs. The 
probability distribution p(n+,nS) develops two maxima at (« + ,«_) = («,0) and (0,«), see 
Fig. Obi), which correspond to fast motion into the plus and minus direction, respectively. 
The cargo stochastically switches between these two maxima and thereby between steady 
plus and minus motion, see Fig.[3lb2). The emergence of these two maxima can be under- 
stood in the following way: When, e.g. because of a stochastic fluctuation, more plus than 
minus motors are bound to the filament (n + > n_), every plus motor experiences the load 
force F c /n + , while every minus motor experiences the larger load force F c /n-. Since the 
unbinding rate increases exponentially with increasing load force according to equation 10, 
minus motors are more likely to unbind from the filament than plus motors. After the un- 
binding of a minus motor, the remaining minus motors experience an even larger load force 
and are even more likely to unbind. As a consequence, an unbinding cascade of minus mo- 
tors happens until no minus motor remains bound. A prerequisite for this unbinding cascade 
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Fig. 4 The motility diagram for the symmetric tug-of-war of N = JV+ = AL = 5 plus and minus motors 
displays the dependence of the cargo motility on the motor force ratio f = F s /Fj and the filament desorption 
constant K = Eq/tCq. The lines separate regions with different locations , n_ ) or numbers of maxima. The 
colors separate the different motility states of no motion (0) (green), fast bidirectional motion ( — h) (yellow) 
and fast bidirectional motion with pauses (red). When the maximum is at (0,0), the cargo is unbound (grey). 
Parameters are for kinesin-1 as in Table [JJ except for the stall force F s and unbinding rate £o, which are 
varied. In addition, the superstall velocity amplitude is scaled as vg = lnm/s F s /pN in order to keep the 
backward slope Vb/F s of the superstall force- velocity equation (5) constant. The crosses labeled (a), (b) and 
(c) correspond to the parameter sets of Fig. [3] 



is that the motors can exert a sufficiently large force to pull off opposing motors from the 
filament, i.e. the stall force F s has to be comparable to or larger than the detachment force 
Fj. For small force ratios / = F s /Fd, the pulling force has only a small effect on motor un- 
binding, so that no instability occurs and the cargo exhibits the blocked motility state (0). 
For large motor force ratio, the transient predominance of one motor type is amplified and 
stabilized by the described dynamic instability, and the cargo spends most of the time in 
states with only one motor type bound. 

(—0+) Fast bidirectional motion with pauses: For intermediate force ratios /, the 
motor number probability p{n + ,nS) displays three maxima, see Fig.[3lcl), a symmetric one 
with n + = n_ corresponding to no motion as in (0), one with n_ = and one with n + = 0, 
corresponding to fast plus and minus motion as in ( — h). The cargo trajectory therefore 
exhibits bidirectional motion interrupted by pauses, see Fig.[3jc2). 

Fig.|4]shows the classification of motility states for the tug-of-war of N = N + = N- = 5 
kinesin-like motors with respect to the relevant motor parameters / = F s /Fd and K = £q/hq. 
As discussed above, the cargo is in one of the three motility states (0), ( — h) or (— 0+). 
The lines in the motility diagram of Fig. [4] separate regions in which the maxima of the 
motor number probability distribution p(n + , «_ ) are located at different motor number states 
(«+,«_). The colors separate regions with different motility states. The detailed calculation 
procedure is as follows: The single-motor parameters are taken to be equal to the kinesin-1 
values as given in Table [T] except for the plus and minus motor unbinding rates £o and stall 
forces F s . The parameter space (F s , Eq) is then explored systematically. All other parameters, 
i.e. Fj, 7Tq, Vf and vg/F s , are kept constant. Note that rather than the superstall velocity 
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Fig. 5 Motility diagrams for the symmetric tug-of-war of N = N + = N- plus and minus motors. Color code 
and motor parameters as in Fig. [4] (e) is identical to Fig. [4] but only the transition lines that separate different 
motility states are shown. The crosses correspond to the full kinesin- 1 parameter set of Table[T] The transition 
lines for N = 1 can be determined analytically as given in (a). 



amplitude vb, the slope vb/F s of the superstall force- velocity-curve ® is kept constantQ 
For each point (F s , Eq), the maxima of the motor number probability p(« + ,n_) is calculated 
as the eigenvector with zero eigenvalue of the transition matrix of the Master equation |Q~||. 
When the maxima between two scanned points change, we zoom in between these points in 
order to determine the transition point more accurately. The lines shown in Fig.|4]consist of 
these points. 

For large desorption constants K, the motors have a low affinity to the filament; therefore 
the number of bound motors at the maxima of /?(«+ , h_ ) in Fig.|4]is low for high K. For very 
high desorption constants K, the maximum is at (n,n) = (0,0), and the cargo is 'unbound'. 
For small force ratios /, the probability distribution p(n + , w_ ) has a single maximum at a no- 
motion state with an equal number of plus and minus motors bound at (n,n) with 1 < n < 5, 
and the cargo is in the no-motion motility state (0) (green). The two neardiagonal maxima at 
(1,2) and (2,1) are also counted as a single diagonal maximum, which in a continuous state 
space would be at (x,x) with 1 < x < 2. For large force ratios /, the motors can generate 
forces large enough to rip off opposing motors since the stall force is large compared to 
the detachment force. This leads to the unbinding cascade described above, and the motor 
number probability has two maxima, one at a state (n, 0) with only active plus and one at 
a state (0, n) with with only active minus motors, with 1 < n < N. In the latter situation, 
the cargo is in the ( — h) motility state (yellow). For intermediate values of /, both types of 
maxima coexist, and the cargo is in the (—0+) motility state with three maxima of the motor 
number probability p(ti + ,nJ) (red). 

Fig.[3]shows the motility diagrams for varying motor numbers N = N + = N- , where only 
the transition lines between the different motility states are shown. If the cargo is carried 
by only 7V + = 1 plus and AL = 1 minus motor, the motility diagram can be determined 



The need for emphasis of this point was kindly pointed out to us by Harry W. Schroeder III of the 
Goldman lab. 



12 



(-+) (0+) (+) 




Fig. 6 Motility states for the kinesin-dynein tug-of-war with parameters as given in Table [TJ (a) A cargo 
which is carried by N + = 5 kinesins and N- = 4 'strong' dyneins is in the ( — h) motility state. The motor 
number probability p(« + ,n_) shown in (al) exhibits two maxima at (5,0) and (0,4), corresponding to fast 
plus and fast minus motion, respectively. As the plus motion maximum at (5,0) is higher than the minus 
motion maximum at (0,4), plus motion is more probable, and the trajectory in (a2) is biased towards the plus 
direction, (b) A cargo which is carried by JV + = 4 kinesins and N- = 5 'weak' dyneins is in the (0+) motility 
state. The motor number probability p(n + ,n_) shown in (bl) has two maxima, one at the fast plus motion 
state (4,0) and one at the no-motion state (1,4). The trajectory in (b2) exhibits plus motion with pauses, (c) 
A cargo which is carried by JV + = 5 kinesins and JV_ = 4 'weak' dyneins is in the (+) motility state. The 
motor number probability p(n + ,n_) shown in (cl) has one maximum at (5,0), which corresponds to fast 
plus motion, as can be seen in the trajectory in (c2). 



analytically. The motor number probability p(n + ,n ) has its maxima either at (1,1) for 
K < exp(-/) (motility state (0)), at (0,0) for K > 1 (unbound), or at (1,0), (0, 1) otherwise 
(motility state ( — h)). The region with the three-maxima motility state (—0+) only appears 
when the number TV = N+ = N- of motors on the cargo is larger or equal to 2, and increases 
for larger motor numbers, see Fig. [5] 

The crosses in the motility diagrams of Fig. [5] correspond to the full set of kinesin-1 
parameters as given in Table[TJ As other molecular motors have parameters of a similar order 
of magnitude, all three motility states (0), ( — h) and (—0+) are within biological parameter 
range. Therefore, the cell can use tuning of the motor parameters in order to drastically 
change the motion of the cargo. 

The locations of the transition lines in the motility diagrams of Figs. Hand [5] are deter- 
mined by the balance of binding and unbinding of single motors under the load force gen- 
erated by the opposing motors. As all unbinding rates are proportional to £o, and all binding 
rates are proportional to tiq, this means that the transition lines should scale with K = Eq/tzq. 
Indeed, when the desorption constant K is scaled by the motor number N, the diagrams of 
Fig lHb)-(f) almost overlay. As a specific example, cargo unbinding occurs when the rate of 
unbinding of the last bound motor, £o, becomes larger than the rate for binding of one motor 
in the unbound state, Ntzq, i.e. when K = Eq/tzq > N. This is indeed the case in Fig. [5] for 
all N. For the other transition lines, this type of reasoning is not exact but leads to a good 
approximation as we will be shown in Section[4] 
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Fig. 7 Motility diagrams for the asymmetric tug-of-war of iV + kinesins and N- (a) 'strong' and (b) 'weak' 
dyneins, with parameters as given in Table [T] The colors indicate the motility states of fast plus motion (+) 
(dark blue), fast minus motion (— ) (light blue), fast bidirectional motion ( — h) (yellow), fast bidirectional 
motion with pauses (—0+) (red), fast plus motion with pauses (0+) (pink) and no-motion (0) (green). The 
motor number probability and trajectory for special choices of N + and N- are shown in Fig. [6] 



Switching between fast plus and minus motors as in the ( — h) motility state has also 
been observed in a motility assay with only one type of non-processive motors, which move 
to the plus and minus end with equal probability [8]. This is formally equivalent to two equal 
teams of non-processive unidirectional motors, with the addition that motors can interchange 
between being a plus and a minus motor. This situation has been investigated theoretically 
in the framework of two-state ratchet models, and also leads to bidirectional motion with 
bimodal velocity distributions [2]. However, because of the non-processivity of the motors, a 
minimal number of 5 motors is required to produce bidirectional motion in this model, and a 
regime (—0+) with plus and minus motion and pauses is not observed. This is in qualitative 
agreement with our results for unprocessive motors, i.e. for high desorption constants K. 
First, a minimum number of motors N > K is required in order to be in the ( — (-) motility 
state. Second, the region of the (—0+) motility states does not extend into the region of high 
desorption constants K, see Fig. [5] 



3.2 Asymmetric tug-of-war 

Bidirectional cargo transport in vivo is typically dependent on two different motor species 
for plus and minus motion. This plus-minus asymmetry can lead to net transport of the 
cargo in one direction. For example, in the motility state ( — h), the plus motion maximum 
of the motor number probability can be larger than the minus motion maximum, which 
leads to longer plus runs compared to minus runs and to net plus motion of the cargo, see 
Fig. |3a). As cargo motion is no longer symmetric with respect to plus and minus motion, 
seven motility states are now possible, corresponding to the seven different combinations 
(+), (-), (0), (-+), (0+), (-0) and (-0+) of the maxima (+), (-), and (0). In the motility 
state (0), the cargo exhibits almost no motion. The motility states (0+) and (+) correspond 
to fast plus motion with and without pauses, respectively, see Fig.^Jb) and (c). Analogously, 
in the motility states (— ) and (—0), the cargo exhibits fast minus motion without and with 
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Fig. 8 Relabeling of the vertex (n + , n_ ) , its four nearest neighbors, and the transition rates connecting these 
five vertices. 



pauses, respectively. The motility states ( — h) and (—0+) describe fast bidirectional motion 
without and with pauses, respectively. 

Fig. [7] shows the motility states for the tug-of-war of N + kinesins and N- dyneins for 
varying motor numbers N + and N_ . As two different values of the stall force have been 
reported for dynein, see Table Q] we show the tug-of-war of kinesins and 'strong' dyneins 
with 7 pN stall force in Fig. [TJa) and the tug-of-war of kinesins and 'weak' dyneins with 
1.1 pN stall force in Fig.|7Jb). In the tug-of-war of kinesin and strong dynein, the opposing 
motors are of similar 'strength', i.e. have similar stall and detachment forces and similar 
desorption constants. In addition, both motors are 'strong' with large stall to detachment 
force ratios. This is similar to the symmetric tug-of-war with strong motors, so that the 
cargo is mostly in the ( — h) motility state of fast bidirectional motion. Only when one motor 
type is much more abundant than the other, this motor types dominates the motion and the 
cargo exhibits fast motion into the direction of this motor type. 

In the tug-of-war of N + kinesins and AL 'weak' dyneins, the dyneins have a much lower 
stall force than the kinesins. Therefore, the kinesins win the tug-of-war for dynein motor 
numbers smaller or of similar magnitude as the kinesin motor number, and the cargo is in 
the plus motion motility state (+). If the number of dyneins is increased, the dyneins act 
as a brake for the kinesins, and the cargo is in the motility states (0+) of plus motion with 
pauses or in the no-motion motility state (0). If the number of dyneins is increased further, 
the cargo is in the motility states ( — h) or (—0+) of fast plus and fast minus motion without 
and with pauses, respectively. For very large number of dyneins, the cargo is in the minus 
motion motility state (— ). The motility state (—0) of fast minus motion with pauses does not 
appear, but is possible for motors with different parameters [29]. 

Fig.|7]shows that for biological motor numbers of a few motors of each type, the kinesin- 
dynein tug-of-war can lead to fast plus motion, fast minus motion, or fast motion in both 
directions. The type of motion can be regulated by choosing the appropriate number of 
motors on the cargo. Many in vivo cargos, which are often earned by kinesin and dynein, 
have indeed been observed to exhibit uni- or bidirectional motion, and also often exhibit 
pauses. Our results for the kinesin-dynein tug-of-war could be directly tested in in vitro 
experiments. 

In Figs.|5]and[7] the (—0+) motility state of fast bidirectional motion with pauses appears 
only if the number of motors N+ +AL is larger than 3. This can be understood by considering 
the state space network of Fig.[TJb). If the number of motors N + +N- is smaller than 4, there 
is simply not enough room for three maxima. This leads to an interesting prediction on the 
number of motors involved in transport. If there are no other causes for pausing, e.g. physical 
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Fig. 9 Constructing the motility diagram for the symmetric tug-of-war of N + = N- = 5 plus against minus 
motors in the sharp maxima approximation: The overlay of the transition lines for (a) the plus and minus 
motion maxima at (n,0) and (0, n) according to equation )22t and for (b) the zero motion maxima at (n,n) 
according to equation fl23t forms the full motility diagram of Fig.[l0] The cargo is unbound (grey) for K> N, 
see equation (24}. 



obstacles, except the molecular motor tug-of-war, bidirectional cargo transport with pauses 
must involve N+ +N- > 4 motors. Therefore, the observation of pauses leads to a lower 
bound on the motors involved in the cargo transport. 



4 Sharp maxima approximation 

Although our tug-of-war model is rather simple, the motility diagrams of Figs.|5]and|7]have 
a complex structure. The qualitative aspects of these motility diagrams can be reproduced 
within a mean field theory for the dynamics of the average numbers of active plus and minus 
motors, as described in [28]. In particular, the mean field calculation reproduces all motility 
states of the stochastic tug-of-war model, with the transitions between the motility states 
becoming saddle-node and transcritical bifurcations. However, it does not lead to analyti- 
cal expressions for the transition lines, and the quantitative agreement with the results of 
the stochastic calculation is poor. In this section, we discuss a rather simple and intuitive 
approximation that leads to analytical expressions for the transition lines of the motility 
diagram. 

We have characterized the motility states of the cargo by the maxima of the motor 
number probability p(n + ,n_) for n + active plus and n_ active minus motors, because the 
cargo spends most of its time in such extremal states. We now focus on situations in which 
the probability is concentrated at the maxima. Thus, we now assume that the probability 
p(n + ,n ) is non-zero only for the maxima locations and for the nearest neighbors of these 
points. For example, the state («+,«_) with < n + < N+, < «_ < N- has the four nearest 
neighbors shown in Fig. [8] The transition rates between these states are given by the rates for 
(un)binding of one plus or minus motor in equations l|7]l and ((SJl and the analogous equations 
for the minus motors. For notational simplicity, we label the extremal state by and the four 
neighbors by 1,2, 3, and 4, and denote the transition rate connecting state i with state j by 
COij, see Fig. [8] The Master equation then reads 

d 4 

^fl) = £(flflto-/b«w), (16) 
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(5,5) motor strength f (5,0),(0,5) 

Fig. 10 The motility diagram for the symmetric tug-of-war of N = N + = N- = 5 kinesin-like plus and minus 
motors, calculated with the sharp maxima approximation. This diagram is an overlay of the two diagrams 
of Fig. [9] i.e. the lines separating the regions with different locations or numbers of maxima are given by 
equations j22H24t . The colors separate different motility states as in the full numeric motility diagram of 
Fig.f4]to which it is very similar. Parameters are as in Fig. [4] except for the superstall velocity amplitude \'b, 
which is set to zero. 



for the extremal state j = = and 
d 



dt 



P, 



-PiCOiO+PoCOQi 



with i = 1 , . 



(17) 



for the four neighboring states. Note that the latter equation contains the 'sharp maxima 
approximation', since it ignores all neighbors of the states i = 1, 2, 3 and 4 apart from state 
0. The steady state solution of equations ([16} and dl7| l fulfills the detailed balance relation 
(OioPi = PoCOoi- The condition for the point j = to be a maximum is 



P > Pi, i-e. m < co i0 



for ;' = 1 , 



(18) 



This is intuitive, stating that the rates C0/o leading into the maximum state should be larger 
than the corresponding outgoing rates 01o;. In terms of the original transition rates as given 
by equations (f7]l and (8), the inequalities dl <S[ i lead to the conditions 
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(19) 



for the eight transition rates between the point (n + , «_ ) and its four nearest neighbors. When 
the transition rates l[7]l and |8]l are inserted, the inequalities d 1 9b become: 
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Biological motors often have a small backward velocity. For example, kinesin-1 has 
vb = 6nm/s, which is small compared to the forward velocity Vp = 1 fim/s, see Table[TJ In 
the limit of vanishing backward velocity vg = 0, the cargo force becomes 

/ max {n + F J+ ,n_F s _} for both «+,«_ >0 
F c (n + ,n-) = < 3 (21) 

[0 for n + = or n_ = 

and the cargo velocity is equal to the single plus (minus) motor velocity v c = vf+ (or v c = 
Vf_) if n+ > 0, n_ = (w + = 0, «_ > 0), and equal to v c = in all other cases. By using the 
relation d2 lb in ( 120b . we obtain four simple inequalities. If these inequalities are fulfilled for 
a given set of motor parameters and numbers, then the state («+,«_) is a maximum state of 
the steady state motor number probability /?(«+,«_). In this way we can obtain all maxima 
for a given set of motor parameters and thereby construct the motility diagram. 

We now illustrate this for the symmetric tug-of-war of /V = N+ = AL motors with identi- 
cal parameters apart from their forward direction. In the symmetric tug-of-war, the relevant 
dimensionless motor parameters are given in equation l !15l >. With these parameters, one ob- 
tains from the inequalities d20l l: 

(i) fast motion maxima at (n + ,n_) = (w,0) and (0,n) for 

fN-n „ f \ N-n+l 

max -,Ns~ n} ) <K< withn= 1.....N. (22) 

\ «+ 1 / n 

The corresponding transition lines separating regions with different maxima locations 
are shown in Fig.|9]a). 

(ii) a no-motion maximum at («, n) for 

N—n f N—n+1 f 

-eT i <K< e'f withn= 1,...,JV. (23) 

n + 1 n 

The corresponding transition lines separating regions with different maxima locations 
are shown in Fig.|9jb). 

(iii) a maximum at (0,0), corresponding to an unbound cargo, for 

K > N. (24) 

For all other points (n + ,n_) with « + ^ n_ and n + , n_ > 0, the inequalities J20| > have no so- 
lution. This means that all maxima that are obtained within the numerical calculation are also 
found in the sharp maxima approximation, and vice versa. Drawing the transition lines be- 
tween the different maxima regions in the (f,K) -plane leads to the motility diagram shown 
in Fig. [10] The maxima (n,0) and (0,n) with only one motor type bound to the filament 
define the ( — h) regime (yellow), while the maxima («,«) with both types of motors bound 
define the (0) regime (green). The overlap region is the (—0+) regime (red). This sharp- 
maxima motility diagram is in good agreement with the exact motility diagram shown in 
Fig. ID The rather simple analytic equations of the sharp maxima approximation can there- 
fore be used to estimate the transition lines. The sharp maxima approximation tends, how- 
ever, to overestimate the existence region of a maximum, and therefore the (—0+) region, 
because the probability for a given vertex is only compared to the probabilities of its next 
neighbors. 

We now apply the sharp maxima approximation to the kinesin-dynein tug-of-war. We 
fix the motor parameters to the values for kinesin-1 and strong or weak dynein as given in 
TableQ] For N + kinesins and N_ dyneins, a state (n + , w_ ) with < n + < N + and < «_ < 
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Fig. 11 Motility diagrams for the kinesin-dynein tug-of-war, calculated with the sharp maxima approxima- 
tion. The motor parameters are for kinesin-1 and (a) strong dynein or (b) weak dynein as given in Table [T] 
except for the superstall velocity amplitude vg, which is set to zero. The motility diagrams from the sharp 
maxima approximation agree well with the exact motility diagrams shown in Fig. [7] 



N- is counted as a maximum of the motor number probability if all four inequalities i !20t are 
fulfilled. With this method, the number and locations of the maxima are determined, which 
defines the motility state for these values of N + and N_ . When N + and N- are varied, the 
motility diagrams shown in Fig. [TT] are obtained. They agree well with the exact motility 
diagrams shown in Fig. [JJ Again, the sharp maxima approximation tends to overestimate 
the existence region of a maximum, and therefore especially the (—0+) region. 

The sharp maxima approximation is based on the simple idea that the probability is 
highly concentrated around single cargo states (n + ,n_). As described in this section, we 
have taken into account only the four nearest neighbors of the center states. It is possible 
to systematically expand this scheme to the 13 next-nearest neighbors etc. However, even 
including only the eight diagonal neighbors (« + ± l,n_ ± 1) already leads to analytically 
intractable expressions since the state network then contains cycles. In addition, this more 
elaborate approximation scheme lead only to relatively small improvements compared to 
the nearest-neighbor approximation, see [28]. Therefore, we focused here on this latter ap- 
proximation, which quantitatively describes the transition lines of the motility diagrams by 
four simple inequalities (20\ . compare Fig.[l0]and Fig.|4]as well as Fig.QT]and Fig. [7] These 
inequalities can thus be used to estimate the motility state of a given experimental system 
without having to resort to numerical calculations or simulations. Furthermore the inequali- 
ties ( |19t show that a transition between motility states occurs when the rates for binding and 
unbinding of a single motor in a maximum state (n+,n_) become comparable. 



5 Summary and discussion 

We have studied a stochastic tug-of-war model for bidirectional cargo transport by two an- 
tagonistic teams of molecular motors. In this model, the motors act independently and are 
coupled only via the mechanical interaction with their common cargo. 

A particularly instructive case is the symmetric tug-of-war of plus and minus motors 
with the same parameters except their forward direction. In this symmetric case of 'equally 
strong' opposing motors, the cargo can be in one of three 'motility states', depending on 
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the values of the motor parameters: For weak motors with low stall to detachment force 
ratio, the cargo is in the no-motion motility state (0) and exhibits only slow motion. For 
strong motors with high stall to detachment force ratio, on the other hand, the cargo is in the 
motility state ( — h) characterized by switching between fast plus and fast minus motion. In 
the intermediate case, it exhibits the (—0+) motility state of fast bidirectional motion with 
interspersed pauses. The latter two motility states correspond to experimentally observed 
bidirectional cargo motion and have been previously attributed to a coordination mechanism 
rather than to a tug-of-war. The fast bidirectional motion is obtained in our tug-of-war model 
via a dynamic instability arising from the nonlinear force-dependence of the single-motor 
unbinding rate. This instability leads to unbinding cascades of one type of motors, so that 
there is a high probability of having only the opposing motor type active at a given time. 

In the asymmetric tug-of-war, e.g. between kinesins and dyneins, the seven motility 
states (+), (— ), (0), ( — h), (0+), (—0), and (—0+) are possible, which correspond to all 
possible combinations of fast plus motion (+), fast minus motion (— ) and pausing (0). 

We have characterized the parameter dependence of the cargo motility in terms of the 
'motility diagrams', which, similar to phase diagrams, describe how the different motility 
states depend on the system parameters. For the symmetric tug-of-war, the relevant param- 
eters are the stall to detachment force ratio /, and the desorption constant K. The transition 
lines separating the motility states (0), ( — h) and (—0+) lie in biologically relevant parame- 
ter ranges of both / and K. This means that the cell can use fine-tuning of these parameters 
in order to achieve large changes of the cargo motility, i.e. in order to regulate its cargo traf- 
fic. For the kinesin-dynein tug-of-war, we have varied the numbers of kinesins and dyneins. 
For biologically relevant numbers of a few kinesins and dyneins, the cargo mostly exhibits 
fast plus motion, fast minus motion, or fast bidirectional motion. However, even in the latter 
case, the cargo may cover a net travel distance, e.g. because runs in the plus direction may 
be longer and/or faster than runs in the minus direction. The net direction and the speed of 
cargo motion can therefore be regulated in two ways: by choosing the appropriate number 
of motors on the cargo, and by tuning the motor properties. 

The motility state (—0+) of fast bidirectional motion with interspersed pauses only ap- 
pears when the total number of plus and motors, N + +AL, becomes larger than 3. If there 
are no other reasons for pausing, e.g. physical obstacles, this leads to the prediction that 
bidirectional cargo transport which exhibits pauses is carried out by at least 4 motors. 

The transition lines separating the different motility regions in the motility diagrams 
have complex shapes, which are hard to understand at first sight. We have used a 'sharp 
maxima' approximation to obtain a simple and intuitive derivation of these transition lines. 
This approximation explains the transitions between the motility states as determined by 
the balance of the motor binding/unbinding dynamics under force. The sharp maxima ap- 
proximation reproduces the complex tug-of-war motility diagrams quantitatively by using 
only four simple inequalities as given by i20\ . It can therefore be used to determine the 
motility state of a given experimental system without having to do simulations or numerical 
calculations. 

The understanding gained from our model allows us to speculate on constraints imposed 
on the numbers and properties of the motors involved in cellular bidirectional transport. In 
order to obtain fast bidirectional cargo motion, the motors must produce large forces - com- 
pared to their detachment forces - in order to initiate the unbinding cascade which leads to 
this fast motion. This dynamic instability is necessary to achieve effective motor coopera- 
tion without the need for an additional coordination complex. This may be one reason for 
the counterintuitive property of kinesin-1 to be able to produce high forces, corresponding 
to a high stall force, without being able to sustain them for a long time, as reflected in the 
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Fig. 12 Lower left comer of the tug-of-war states network of Fig. |TJb), with the addition of state U repre- 
senting an unbound freely diffusing cargo. 

lower detachment force. The dynamic instability also provides an explanation for kinesin's 
intermediate processivity of 'only' about a second before unbinding: longer unbinding times 
would slow down the unbinding cascade and lead to lower velocities and longer pauses. In 
order for cargo motion to be easily controllable, a small change in the motor parameters 
should lead to a large change in cargo motility. Such a response behavior is obtained in 
our model when the motor binding and unbinding rates are of similar order of magnitude 
corresponding to a desorption constant of the order one. 

Finally, why should the cell use two teams of opposing motors on a single cargo at 
all? Using only unidirectional motors on a single cargo is not sufficient to maintain cellular 
transport in an efficient way. As the cellular microtubule cytoskeleton is typically isopolar, 
plus end cargos would accumulate at the microtubule plus ends, and minus end cargos at 
the microtubule minus ends. Even if motors are interchanged at the microtubule ends, which 
would solve the cargo jamming problem, it would still leave the problem how the plus (mi- 
nus) motors reach the minus (plus) ends in order to start their journeys. Having motors of 
opposite directionality on the same cargo is an elegant solution to this problem. In addition, 
bidirectional cargos can easily be regulated. By using the sensitivity of the cargo motility 
to the motor properties, the cell could influence the cargo motion easily by fine-tuning the 
parameters of the two motor species. 
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A Cargo unbinding 

In this appendix, we show that the bound motion of the cargo is independent of the rates for cargo unbinding 
and binding. 

We have characterized the cargo state by the numbers n + of plus and ;i_ of minus motors that link it 
to the filament. Cargo unbinding can only occur when no motors link the cargo to the filament, i.e. from the 
state n + = re_ = 0. This state can be split into two states, see Fig. 1121 an unbound diffusive state U where the 
cargo is far away from the filament, and a state (0,0) where the cargo is close to the filament as if bound by 
motors (but no motor is bound). The cargo is in the (0.0) state, if for example the last motor linking it to the 
filament has just unbound a short time ago. From state (0,0), the cargo can 'rebind' to the filament because 
of binding of a plus or a minus motor. Alternatively, the cargo can 'unbind' from the filament by diffusing 
away into the unbound state U. This is described by the unbinding rate eu. The binding rate Ku from state 
U to state describes the probability that the freely diffusing cargo comes close enough to the filament for a 
motor to bind. The rates 7C\j and £\j relate to diffusion of the cargo in the surrounding solution and therefore 



21 



depend on the geometry of the system, e.g. the location of the filaments or the volume and viscosity of the 
surrounding medium. 

The rates leading from the (0,0) state to the states (1,0) and (0, 1) with one bound plus resp. one bound 
minus motor are the rates for the binding of the respective motor when the cargo is close to the filament. These 
are the rates 7T + (0,0) = NlZo+ resp. 7T_ (0,0) = Nhq- used in Section[2] Actually, all the rates X+i-) («+,»- ) 
use the single motor rates Hq± which are the rates for the binding of a plus resp. minus motor to the filament 
when the cargo is close to the filament. 

Let P(n + ,n_) denote the steady state probability of n + plus and n_ minus motors linking the cargo to 
the filament, and P(U) the steady state probability of being in the diffusive unbound state U. As experiments 
only monitor bound cargos, we consider conditional probabilities, conditioned on the cargo being bound: 

p(n+,n-) = P(n + ,n-)/[l-p(U)} (25) 

for < n + < N + and < n_ < N-. The state (0,0) is included here, since in this state the cargo has not 
diffused away yet and appears as bound to an experimenter. 

We now show that in the steady state, these conditional probabilities do not depend on the rates 7Tu and 
£u . In order to do this we consider the calculation of the steady state via the diagram or Kirchhoff method. This 
method is similar to the Kirchhoff rules for electrical networks and is reviewed in [33]. First, one constructs 
the complete set of partial diagrams for the state space network. A partial diagram is a diagram with the 
maximum number of lines that can be included without forming a closed cycle. From these partial diagrams, 
one constructs the directional diagrams for each state n of the network: add arrows to each partial diagram 
flowing towards state n. In our case, n can be one of the states (n + , «_ ) with < n + < N + and < «_ < JV_ , 
or U. Each directional diagram represents the product of the rates of the arrows of the directional diagram. 
Let II" denote the product of rates associated with the ith directional diagram of state n. Then the steady state 
probability of state n is given by 

PW = E"nf/s. (26) 

Here, S = Y,„ £; n" is the sum of all directional diagrams. denotes the sum over all directional diagrams 
i of state n, and the sum over all states. Because of the equality 

l-P(U) = (S-l,'/n^)/S, (27) 

the conditional probability to be in state n under the condition not to be in state U is 

Pin) Y"n" 

pin) = { - = u ' „ for n + U. (28) 

Now consider the directional diagrams II" in more detail. They all contain the line U — (0,0), as this line is 
never part of a cycle. Any directional diagram Yl^ 1 for state U contains the rate £[/, but not the rate 7lu, as all 
paths in these directional diagrams must end in state U. On the other hand, all directional diagrams J7" with 
n y^U contain the rate Tlu but not £y (otherwise there would be a path ending in U, but all paths must end 
in n). Thus in equation j28l . £u does not appear (the terms in the total sum 5 in the denominator containing 
Ea are substracted) and the rates Tlu appear exactly once in every term in the numerator and denominator and 
thus cancel out. Therefore, p(n) does not depend on the rates Ky and £[/. 
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